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The Dirichlet Problem for the fractional 
p-Laplacian evolution equation 

by Juan Luis Vazquez * 


Abstract 

We consider a model of fractional diffusion involving the natural nonlocal version 
of the p-Laplacian operator. We study the Dirichlet problem posed in a bounded 
domain U of with zero data outside of U, for which the existence and uniqueness 
of strong nonnegative solutions is proved, and a number of quantitative properties 
are established. A main objective is proving the existence of a special separate 
variable solution U(x,t) = F(x), called the friendly giant, which produces 

a universal upper bound and explains the large-time behaviour of all nontrivial 
nonnegative solutions in a sharp way. Moreover, the spatial profile F of this solution 
solves an interesting nonlocal elliptic problem. We also prove everywhere positivity 
of nonnegative solutions with any nontrivial data, a property that separates this 
equation from the standard p-Laplacian equation. 
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1 Introduction 


In this paper we consider a model of fractional diffusion involving nonlocal operators of 
the p-Laplacian type, and we describe the main properties of the solutions of the evolution 
process posed in a bounded domain hi of ~R N with zero Dirichlet conditions outside of fi. 

Our main contributions are the following: after posing the problem in the frame¬ 
work of semigroups generated by maximal monotone operators, showing existence and 
uniqueness of so-called strong solutions, and establishing the basic quantitative proper¬ 
ties, we prove the existence of a special separate-variable solution U(x, t ) of the form 
U(x,t ) = F(x), that we call the Friendly Giant. Special solutions have played a 

great role in the development of very many areas in nonlinear analysis and applied mathe¬ 
matics. This happens in particular in the area of diffusion and heat transport, where there 
are some famous instances like the Gaussian kernel, the Barenblatt solutions, different 
traveling wave solutions, and other interesting examples with a lesser influence. 

Our special solution plays a similarly important role for the problem we are addressing: 
it produces a universal upper bound for the whole class of nonnegative solutions of the 
homogeneous Dirichlet problem. This means that the whole class of solutions is uniformly 
bounded at any given time t > 0, and the universal bound decays to zero as time grows. 
Actually, the accuracy improves as t —» oo : U(x,t) gives the first approximation the 
asymptotic behaviour of the all nonnegative solutions. Finally, we show that the spatial 
profile F(x) of this solution solves an interesting nonlocal elliptic problem. 

Let us be more precise: our aim is to study the equation 

».(*■«)+/ a-!) 

Jrn \x - y\ N +*P 

where x G C M. N , N > 1, $(^) = c\z\ p ~ 2 z, 1 < p < oo, and 0 < s < 1. The constant 
c > 0 is unimportant. The integral operator may be called the (p, s)-Laplacian operator, 
or more precisely the s -fractional p-Laplacian, our notation is C P)S . In this paper we 
cover the range 2 < p < oo. Note that for p — 2 we obtain the standard s-Laplacian heat 
equation, u t + (—A ) s u = 0, about which much is known; on the other hand, it is proved 
that in the limit s —> 1 with p ^ 2, we get the well-known p-Laplacian evolution equation 
dtu = A p (u), after inserting a normalizing constant. Due to this, the notation (—A)* is 
used in some papers for C P)S . 

We recall that the standard p-Laplacian with p > 2 is a well-known example of nonlinear 
diffusion of degenerate type, it is important in the study of free boundaries and as a 
benchmark for more complex nonlinear models, cf. the monograph [5X1140], The fractional 
version we are going to study has attracted attention only recently with the current focus 
on the analysis of nonlocal diffusion models of fractional Laplacian type, see the surveys 
[13156]. References to previous work on equation (II.Ill and its stationary versions are 
[3 EE [27], [39] . There is much current activity on this topic, see [26] [35] [38] El US [46] . 

We will consider the equation in a bounded domain D C with initial data 

u(x, 0) = u 0 (x), x G D, (1.2) 
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where Uq is a nonnegative and integrable function. Moreover, we impose the homogeneous 
Dirichlet boundary condition that in the fractional Laplacian setting takes the form 

u(x, t) = 0 for all x E ~R N , x qL and all t > 0. (1.3) 

When we apply operator C P}S on the set of functions that vanish outside of we may 
be more precise and use the notation C P)S ,th but this longer notation will not be used 
since there is no fear of confusion. We call Problem DP the problem of solving (II.ip with 
conditions (II.2p and (II.dp . More precisely, we should write DP(£ ps ;Q). 

Outline of the paper 

In Section [2] we gather the necessary preliminaries and notations, we analyze the differ¬ 
ential operator C PtS , and we state the basic existence and uniqueness results that follow 
from standard theories once the problem is correctly posed. We also derive the basic 
properties of the evolution semigroup corresponding to Problem DP. 

In Section |3] we state and prove some of our main results: the existence of the special 
separate variable solution U(x, t ) that we call the Friendly Giant, and its role as universal 
upper bound, Theorems 13. II and 13.21 A historical source reference for this type of result is 
the paper by Aronson-Peletier [3] on the Porous Medium Equation, 1981. The work was 
continued in where the curious name originated. A complete account of the results, 
for positive and sign-changing solutions was given in [53] . see also [ 55] . 

The special profile F(x) of our solution U(x,t) is the unique positive solution (ground- 
state) of an interesting nonlocal elliptic problem, see formula (13. 2 p in Theorem 13.11 

Section [H deals with the asymptotic behaviour. We establish the correct rate u(x,t ) = 
0(f -1 /( p ” 2 )) and the sharp form u(x, t) —* F{x) in Theorem 14.11 Consequently, 

we prove that the asymptotic behaviour does not depend at all on the solution (it is 
universal), as long as u$ is nonnegative, nontrivial and integrable. 

In Section [5] we investigate whether the solutions have the Unite propagation property 
of the standard p-Laplacian or, on the contrary, the nonlocal effect produces positivity 
at all points and for all times. The latter option is concluded: long distance effects win 
and nonnegative solutions are positive everywhere. 

Finally, in Section [h] we give an improved version of the asymptotic convergence theorem, 
with sharp rate of convergence, under a boundary regularity assumption similar to the 
one established in [47, 33] for p = 2. We state this assumption as an important conjecture. 

The paper concludes with a long list of comments and open problems. 

Fractional Laplacian operators, motivation and related works 

Let us recall that fractional Laplacian operators and related nonlocal operators of dif¬ 
ferent types have attracted considerable attention in recent times because of the interest 
in the applications and also because of the rich mathematical theory that has arisen. 
Indeed, interest was fueled by a number of applications in various fields like continuum 
mechanics, stochastic processes of Levy type, phase transitions, population dynamics, 
image processing, finance, and so on. The original models were mostly linear, indeed 
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the linear fractional heat equation has a long tradition in Probability since fractional 
Laplacians are the infinitesimal generators of stable Levy processes mum- However, 
a decade ago intense work began on nonlinear models, featuring both stationary and 
evolution models. Among the first, we mention the work on nonlocal obstacle models, 
m Q23 , on semilinear elliptic problems of nonlocal type [HJ [48] j49], on phase transitions 
[5Tj . and the work on nonlocal perimeters PH ESQ. 

Evolution problems of nonlocal diffusion type originated from material science and sta¬ 
tistical mechanics, and led to much work on fractional porous medium models, see for 
instance mmmmmmummmm, the work was summarized in 'oTi.TKI. The 
Dirichlet Problem for a version of the fractional Porous Medium in a bounded domain has 
been recently treated by the author in collaboration with Bonforte and Sire, [9] [IT]. The 
close parallel that has been found between the theories of the standard Porous Medium 
Equation and the p-Laplacian equation, accredited in publications like [31 » Ml SI], and 
the progress achieved in the study of fractional porous medium models makes it natural 
to investigate the fractional p-Laplacian equation (11.11) . 

We will be inspired in the techniques developed to treat the Dirichlet Problem for 
the standard Porous Medium Equation, as explained in detail in [S3], where references 
to related work are found. The main ideas involved are rescaling, existence of special 
solutions, a priori estimates, and monotonicity arguments. The rescaled orbits converge 
to stationary states which solve a nonlinear elliptic problem. Such general techniques 
were then applied to the evolution p-Laplacian equation, and then the fractional Porous 
Medium. 

The existence of a certain special solution is crucial in our analysis. Indeed, almost no 
special solution is known for equation (11.11) . In this respect we should mention Lindgren- 
Lindqvist’s [32j who End the first eigenvalue and eigenfunction of the operator £ p s on a 
bounded domain and describe its properties, in particular the limit as p —* oo. 

Finally, the unpublished paper [4Tj contains the derivation of the basic existence and 
uniqueness results for different initial and boundary value problems in the language of 
semi-group theory. Together with [36 , 39] , it allows us to shorten a bit our preliminaries 
section. In fact, the basic theory falls quite well into the scope of classical nonlinear 
semigroup theories for monotone or accretive operators, |T311 25], 

Reminder and notations. In the sequel p > 2 and 0 < s < 1 will be fixed. Only 
nonnegative solutions will be considered. The space-time domain of the solutions is 
denoted by Q = x (0, oo). Also, R denotes the closure of 0, and d(x) = d(x,dQ) 
denotes the distance from a point x G R to the boundary dVt. The notation f + denotes 
the positive part of a function, i. e., max{/(x),0}. We write both d t u and u t for the 
partial derivative in time, either for clarity or brevity. 
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2 Preliminaries 


• Given numbers p G (1, oo) and s G (0,1) we define the Gagliardo functional for smooth 
integrable functions defined in as 




1 

p 


u(x) — u(y) \ p 
\x - y\ N+s P 


dxdy. 


( 2 . 1 ) 


By extension it defines a norm in the fractional Sobolev space (also called Gagliardo or 
Slobodetski space) Vh s,p (M Ar ). Indeed, the norm is given by 


ru 


W S ’P 


= Jp,s(u) + / \u(x)\ p dx. 


( 2 . 2 ) 


We refer to |32j for information about these spaces and their properties. It is known that 
as s —> 1 the space W s,p becomes W 1)P , see [12]. 

Functional J P}S is a convex, lower semi-continuous and proper funcional, and it has an 
associated Euler-Lagrange operator, given in its weak form by the expression 



u(y)\ p 2 (u(x) - u(y)) {<f>(x) - <t>{y)) 
\x — y\ N+s P 


dxdy 


(2.3) 


where </> is any smooth variation in IF p,s (M Ar ). This defines the operator C p>s as 


£ p , s (u)(x) = 2 p.v. 


u(x) ~ u(y)\ p 2 (u(x) - u(y )) 

\ X _y\N+S P a y ’ 


(2.4) 


where “p.v.” stands for the principal value of the integral. Thus defined, it is a positive 
operator, that corresponds for p = 2 to the standard definition of (— A) s (up to a constant 
that we do not take into account). Caveat: we can find in the literature the definition 
with the opposite sign, to look like the usual Laplacian, but we will stick here to this 
definition as a positive operator, following the tradition in fractional Laplacian operators. 

• When restricted to the Hilbert space L 2 (R N ), jC PtS becomes a maximal monotone oper¬ 
ator, as the sub-differential of functional ,J p , s . It has a dense domain 

D(C p<a ) = {ue L 2 (R n ) : J ps {u) < oo, £ P}S u G L 2 (R N )} C L 2 (R N ). 


This means a number of properties for the evolution equation d t u + C P)S u = 0. According 
to the general theory in monograph [T3], the usual mild solutions of the evolution equation 
provided by the semigroup theory are indeed strong solutions, which means that for every 
u 0 G L 2 (M^) the evolution orbit u{t) : C([0, oo) : L 2 (R N )) is differentiable for a.e. t > 0 
and du/dt = —C P)S u G L 2 (R N ). Besides, 


C 

N ) < ~ ||«o||z,2( K 2V). 


(2.5) 


4 





Moreover, a number of properties of £ ps are similar to the standard p-Laplacian operator, 
— A p . The convergence as s —> 1 of the solutions of C p ^ s = f to the solutions of —A p u = f 
is proved in . As a token of that the similarity, it is shown that £ ps is contractive in all 
L 9 -norms, 1 < q < oo, and this makes it into an m-completely accretive operator in the 
sense of |4j, see proofs in [41] . Consequently, the map St : uq h > u(t ) is an L 9 -contraction 
in L q (R N ) for every 1 < q < oo (and every t > 0). We can also use the Implicit Time 
Discretization Scheme to generate the evolution solution by solving the elliptic problems 

h £ p, s u T ii f 

with h > 0 a small constant (the time step), and applying the Crandall-Liggett iterative 
scheme [29]. 

• Let us now consider the problem posed on a bounded domain 12 C M A . For simplicity we 
will assume that 12 has a smooth boundary but this is condition can be relaxed. In order 
to take into account the homogeneous Dirichlet boundary conditions in the fractional 
setting, we restrict the operator to the set of functions 1F 0 S ’ P (12), defined as the closure 
of C“(12)0 with respect to the norm induced by J P:S according to (12.2D . The formula for 
the sub-differential is the same, only the functional space changes. Indeed, formula (12. 4 p 
still applies as the defining formula with x G 12, but integration in y is performed all over 
M. N . This means that we can feel the influence of the value u(y) = 0 imposed for y ^ 12 
on the diffusion process inside 12, since the formula involves the differences u(x) — u ( y ). 

We thus finally obtain a maximal monotone operator defined in a dense domain in L 2 (12) 
that is also m-completely accretive, see more in [4T]. We should denote this bounded 
domain operator as £ PtS ,n to avoid confusions with the operator defined for all x G M A? , 
but we feel that this extra care is not needed, so jC P)S will be used in the sequel when we 
work on a bounded domain. In any case, there are many ways in which the difference 
between the action of C p>s in the whole space and on a bounded domain is reflected, a 
very clear case is made in the first comment of Section [7] concerning the lack of mass 
conservation for problem DP, for which we obtain a simple formula. 

• We get the following existence and uniqueness result. 

Theorem 2.1 For every w 0 G L 2 (f2) there is a unique strong solution of Problem DP 
such that (i) u G C([0,oo) : L 2 (f2)) ; (ii) for every t > 0 C P)S u(t ) G L 2 (f2), and (in) 
equation (II.ip is satisfied for a.e. x G 12 for every t > 0. Moreover, for every two 
solutions U\,U 2 with initial data woi,wo 2 we have 

II (ui(t) - li 2 (t))+||l < ||(«01 - « 02 ) + ||l , ||«1 ( t ) - U 2 (t)\\ q < ||w 0 i - u 02 \\q (2.6) 

for every 1 < q < oo and t > 0. 

Compare with Theorem 2.5 of |4lj . The theory can be applied to solutions with any sign. 
In the sequel we will deal with nonnegative solutions, u(x,t),uo(x) > 0. We also point 

lr These functions are naturally extended by zero outside of f2. 
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out that some references introduce the concept of viscosity solutions, like [27J ;36| [39], 01] , 
but we will not need this concept here. 

The class of initial data can be extended to L 1 (h2) using the iWcontr activity. The 
resulting class of solutions is the same for t > 0 since we will prove that all solutions are 
uniformly bounded for t > r > 0, hence u(t) G L°°(Q) C L 2 (Q). 

Here are some properties of the solutions that will be used below. 

• Scaling property: a simple inspection shows that C P)S (Af(Bx )) = A p_1 B ps (C P)S f)(Bx) 
for any A, B > 0. When applied to the strong solutions of the evolution problem, it means 
that when u(x,t) is a solution defined in the domain 9, then 

u(x,t) = Au(Bx,Ct) (2.7) 

is a solution on 9' = {x : Bx G 9} on the condition that C = A p ~ 2 B ps . If we want to 
keep the domain fixed, then we need B = 1. Notice that we also have C p%s (—f (x)) = 
— (yjC,p S f) (x), so that — u(x,t ) is also a solution if u is. 

• The collection of nonnegative solutions has the property of almost monotonicity in 
time 

'll 

dtU - - (2 ' 8) 

valid a.e. in 9 for every t > 0. The proof is based on the Benilan-Crandall homogeneity 
argument [2]. Let us recall it briefly: scaling says that if u\(x, t ) is a solution, also 
u\(x,t ) = Xu(x, \ p ~ 2 t) is a solution. But if A > 1 we have ma(^,0) > Ui(x, 0) > 0, hence 
by the maximum principle we get u\(x,t ) > ui(x,t). Differentiate in A for A = 1 to get 
(p — 2 )tu t + u > 0. A related homogeneity argument proves that 

\\dtu\\ q < for evei A Q > 1 • ( 2 -9) 

We ask the reader to complete this proof, see a similar argument in (532 pg. 185]. 

• When fl is a ball, say Br(Q), and uo is radially symmetric, so is the solution u(x,t ) 
as a function of x . This is an immediate consequence of uniqueness and the rotation 
invariance of the operator. 

• Comparison of solutions in different domains applies: if Ui, i — 1,2 are two solutions 
in domains 9* such that fli C O 2 and we also have Ui(x, 0) < U 2 {x, 0) for all x G fli, then 
U\(x,t) < U 2 (x,t) for all 2 : G 0| and t > 0. 

• A usual property of nonlinear evolution processes of diffusion type is the so-called 
smoothing effect, whereby general initial data in a Lebesgue space L p (h2) will produce 
bounded solutions for positive times, with an L°°(Q) that depends on the norm of the 
initial data. In the present case, much more will be proved, the bound will not depend 
on the initial data. This is what we study next. 

• Before going into that, let us state a result hat is not evident at all. 
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Theorem 2.2 Let u be the solution of problem DP corresponding to nontrivial initial 
data Uq E L 1 (Q), u 0 > 0. Then u(x,t ) is positive in the whole domain Q = ft x (0, oo), 
and uniformly positive on compact subsets. 

The proof of positivity will be postponed to Section [5] This result eliminates the possible 
existence of free boundaries, which were a typical feature of the standard p-Laplacian 
equation. For the fractional porous medium equation, a similar infinite propagation 
result holds in the model studied in [121 01], but free boundaries appear in the model 
studied in |22. [23] . 

• The question of regularity of the stationary problems £ PiS (w) = / has been studied 
recently in [20 133] . Theorem 1.1 of the latter paper states that when / E L°°(ft) then 
u E C“(ft) for some a > 0. Due to the explained regularization effect of the evolution 
problem DP we have Ut E L°°(ft) uniformly for every t > r > 0, hence the solutions of 
Problem DP are uniformly Holder continuous w.r.t. the space variable for t > r > 0. C a 
regularity is directly obtained by Sobolev embeddings when sp > N. 


3 The Friendly Giant. Boundedness of the solutions 

We will prove the following two results 

Theorem 3.1 Let p > 2. There exists a special function U(x, t) of the form 

U(x,t) =t~ 1 /^F(x) (3.1) 

that is a strong solution of Problem (II.II) - (II.3D for every t > 0. The profile F is positive 
in ft and vanishes on <9ft. It is a multiple of the unique positive solution of the stationary 
problem: 

f E W^R*), = f { x ) in ft, f(x) = 0 for x E M. N \ ft . (3.2) 

Moreover, F is C a -Hdlder continuous for some a > 0, and it behaves near the boundary 
like 0(d(x) s ). 

Note that F depends on ft, p and s. As a side remark, the special solution (13.ip satisfies 
dtU = —U/((p — 2 )t), which means that the a priori estimate (I2.8|) is indeed optimal 
(with optimal constant). 

This special solution will be called the Friendly Giant as we mentioned above. Strictly 
speaking, U(x, t ) is not an admissible solution of Problem DP since U(x, 0) = +oo for 
every x E ft. The easy remedy is to apply a small time delay and consider the solutions 
U(x,t + h ) for h > 0. On the other, it is precisely the unbounded initial data what makes 
it an interesting object in the theory. Thus, we will prove that it is a universal bound for 
all nonnegative solutions of the Dirichlet problem, as the following result shows. 
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Theorem 3.2 For every nonnegative solution u(x,t) with initial data in some L p {Fl), 
1 < p < oo, we have 


u(x,t) < U(x,t) := t 2 ^F(x). 


(3.3) 


3.1 New variables 

In order to prove the theorems, we start by making the change of variables 

v(x, r) = (a + tY^-^ufx, t ), r — -log(a + t) (3.4) 

p — 2 

with a constant a > 0. Making the calculations and using the scaling property of £ PiS , 
the new equation is then 

d T v + (p - 2)£ PtS (v) = v . (3.5) 

This is called the rescaled equation. If a = 0 we have a small problem with the initial data 
since t — 0 implies r = —oo, so that the solution v(x, r) is an eternal solution, defined 
for — oo < r < oo. But on the other hand, it has a clear advantage: the property of time 
monotonicity (12.8p is equivalent to saying that d T v > 0 in the sense of distributions (at 
least). This is an important tool at our disposal. 

3.2 A uniform bound 

Lemma 3.3 There is a bounded function W(x) > 0 such that for every solution of 
Problem DP we have 

u(x,t) <t~ 1 ^ p ~ 2 ' ) W(x). (3.6) 

Proof. It is a question of finding a suitable candidate. Here the candidate W(x) is 
assumed to have a top mesa of height 1 (i.e., a maximum value W(x) = 1) in the ball 
of radius 2 R (which is going to be large), W is radially symmetric and decreasing for 
|x| > 2 R, and decreases to zero as |a;| -4 oo in a way to be chosen. The use of the 
integral formula (12.4j) shows that <3>(:r) = £ P)S W(x) is positive and smooth for |x| < 2 R 
(i.e., at all points of maximum), and it is smooth outside. Hence, for \x\ < R we have 
£p,sW{x) > C > 0. Since IT < 1 we get 

C P)S W{x) >CW{x) for |x| < R. 

Putting XW instead of W we can fix the constant C to any given value, say C — 1. 

Remark. We do not know if £ PtS W(x) is positive for |x| > 2R, it depends on the tail; 
for compact support it is negative with decay 0(|x|^ Ar+sp ). This does not look so good, 
but such delicate information is not needed in our next argument, that goes as follows: 

• We perform the comparison step. We assume that C Rr( 0). We write 

d T (v - W) + £ p , s (v) - £ p , s ( W) = v- £ P , S (W) for x e D. 
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We multiply this equation by p(v — W) where p G C 1 (M) with p(s) = 0 for all s < 0, 
p(s) = 1 for s > a > 0, and p'(s) > 0 for 0 < s < a. If j is the primitive of p with 
j(0) = 0, we then get 


" 7 “ f j{v — W)dx + 
dr .Ln j v ' 


{£ p ,s(v) - £p, s (W))p(v - W)dx = / (v - £ P)S W)p{v - W)dx. 


Since v(x) = 0 for x ^ fl and W(x) >0 everywhere, we have p(v(x) — W(x)) = 0 for 
all x 0, so that the integrals can be taken indistinctly over or over M w . The second 
term of this display can be written in two parts. The first part is 

/' f (\v(x)-v(y)\ p - 2 (v(x)-v(y)))(p(v(x)-W(x))-p(v(y)~ W(y))) 


x 


y I 


N+sp 


dxdy , 


and similarly the second part using £ p<s (W) instead of £ PtS (v). Subtracting them, we see 
that the whole term is positive. Finally, we replace the last integral by 


(v- £ PtS W)p(v - W)dx < (v- W)p(v - W)dx , 

Jfi J n 

using the pointwise inequality £ PtS W > W in Putting all together, we get 

-y- [ j(v — W)dx< f (v — W)p{v — W)dx . 
dr J r n J n 

We may let p(s) tend the Heaviside function so that j(s) —>■ s + . Putting now X(r) : = 
f(v — W ).|_ dx, the limit of the last inequality gives 


X{t) > 0, — <X(r). 

If the initial data are zero, then the solution is always zero. But this means that v(x, r) < 
W(x) for all x G Vl and all r > 0. This is equivalent to saying that t~ 1 ^ p ^W{x) is an 
upper bound for the corresponding solution of Problem DP. 


• Can we ensure that X(t) = 0 at some given time? It is not difficult when the DP starts 
with a function u 0 is bounded and has compact support in D. We proceed as follows we 
want to compare u(x,t) with (t + £)~ 1 ^ p ~ 2 ' l W{x) before doing the change of variables 
(i.e., the rescaling). The comparison is true at t = 0 if e is small enough since W is 
uniformly positive in D. We now pass to the rescaled formulation with a = e so that 
t — 0 corresponds to some To > — oo. Since W(ro) = 0, we can apply the previous result 
to conclude that X(r) = 0 for all t > Tq- This means that 

u(x, t) < (t + eY l ^ p ^W(x) < W(x), 

which ends the proof. 

• We want to prove that the upper bound holds for any solution of Problem DP with 
data uq G L 1 (h2). In this case we approximate uo by compactly supported and bounded 
functions like in the previous step, get the uniform estimate for the approximations, and 
pass to the limit using the contraction property of the semigroup. The proof is concluded. 

□ 
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3.3 Alternative method using eigenfunctions 


In order to prove the upper bound, we may use as starting barrier the eigenfunction 
results that have been obtained by Lindgren and Lindqvist in p|j. They consider the 
minimization of the fractional Rayleigh quotient 


^p,s(0) ; — 


1 <i>(y) - (f>(x)\ p 

\x — y\ N+s P 


dxdy 


\4>(x)\ p dx 


among all functions <fi in the class 0^0. They assume that p >2 and 0 < s < 1. 

They also write ap instead of N + sp, so when 0 < s < 1 we get N < ap < N + p 
(this change of notation is not important here). They call the inhmum of J PjS (f2) the 
first eigenvalue, Ai. It is positive and depends on p, s and fl Actually, it increases as we 
shrink by a scaling argument. We are interested in Theorem 5, Section 3, of [39], that 
says that: 

Existence of the first eigenfunction. There exists a 11011 -negative minimizer f G 
W s,:P (f2) with f ^ 0 in and / = 0 in R N \Q. It satisfies the Euler-Lagrange equation 

w^i/r 2 / ( 3 . 7 ) 

in the weak sense with test functions 0 G C'^°(f2). If sp > N, the minimizer is in C' 0,/3 (M JV ) 
with (3 = s — N/p. 

• We use the previous result to construct a supersolution of the eigenvalue problem 

C P:S (H) = cH 

inside a ball Br by using the first eigenfunction / of the minimization problem posed in 
a larger domain, like the ball B 2 r. This will be a way to obtain the upper for v with a 
more explicit construction. 

In this way we may retrace our steps and prove the existence of the universal upper 
bound if p is large, more precisely sp > N. [39] does in that case a viscosity theorem and 
proves uniqueness and positivity of the first eigenfunction. But for smaller p we have the 
regularity and positivity results that have been proven recently. 


3.4 The bounded limit profile F{x ) 

Let us go back to our rescaled solution v(x, r). We have monotonicity time of the whole 
orbit for —00 < t < 00 . We also have a uniform a priori bound v(x , r) < W(x) valid for 
all x G and all r G R. Therefore, we are allowed to pass to the limit t —> 00 and get 
the expression 

lim v(x , r) = F(x) (3.8) 

t—too 
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with F(x) < W(x). Since we also have a uniform bound of the form on d T u and C PtS u, it 
is not difficult to prove that the equation (j3.5j) becomes 

£pA F ) = P F , H= —^ , (3.9) 

at least in the weak sense. The factor fi is not essential. It can be eliminated by a simple 
rescaling that uses the different homogeneities of the right-hand and left-hand side; is we 
put F = cF\ and c p_2 = /i, then C PjS (Fi) = F\. 

Finally, since Theorem 12.21 asserts that u(x,t) is positive everywhere, so is v(x,t). F(x) 
is a monotone increasing limit of v, so F(x) > 0 for every x G hi. We can now apply the 
regularity results of [2(Jl 'Toj to conclude that F is C a Holder continuous for some a > 0 
and that it has the stated behaviour at the boundary (see [35| Theorem 4.4]). 

3.5 Universal bound and uniqueness of the profile 

By means of passage to the limit of a particular nontrivial orbit we have obtained a 
Friendly Giant of the form U(x,t ) = F(x) that is an upper bound for the whole 

orbit u(x, t ) with initial data Uq(x, t) from which we started. But the argument of Lemma 
13.31 applies to any other solution u\{x, r) with initial data 'Uoi(x) > 0 using now F as an 
upper barrier in the place of W. In this way we conclude that t~ l ^ p ~ 2 ' > F(x) is a universal 
upper bound for any orbit. 

Following this argument we can prove the uniqueness of the strong nonnegative solution 
of problem (13.51) . Suppose we have a second solution F > 0. We consider the function u 
solution of (II.ip with data F at t — 0. By the previous argument we have 

u(x,t) = (t + Fi(x) < U(x,t) =t~ 1 l (jp ~' 22> F(x) 

passing to the limit we conclude that F\ (x) < F(x). The roles can be reversed, hence 
Fi — F The proof of Theorems 13.11 and 13.21 is done, pending the proof of positivity of the 
solutions. □ 

Remarks. (1) We have proved that equation (13.51) for v(x, r) has just one positive 
stationary solution. However, it has many non-stationary ones, derived from it and 
following the formula 

V (x, r; a) = (1 + -) 1 /(p~ 2 ) F(x) . 
t 

All of them tend obviously to F as t — * oo. 

(2) If we expand the domain by means of a homothetical scaling from H to AH = {Ax : 
x G H}, then we get the following formula for the corresponding profile 

F xa {x) = \^^F a {x/\). (3.10) 
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4 Asymptotic behaviour 


Here, we establish the asymptotic behaviour of all nonnegative solutions of the Dirichlet 
problem. 

Theorem 4.1 For every nonnegative solution u(x,t ) of the DP we have 

lim t 1 /(p _2 )|| M (x, t) — U(x,t)\\ q = 0 . (4.1) 

t—>oo 

for every L q norm with 1 < q < oo. In other words, ^/(p -2 ) u(x,t) tends to F(x) as 
t —y oo uniformly in 12. 

Proof. The statement is equivalent to saying that v(x,t) converges uniformly to F(x). 
Now, the time monotonicity and the boundedness of F allows us to use the monotone 
convergence theorem to conclude that v(-,t) converges to F(-) in L 1 (l2). This and the 
uniform boundedness imply the convergence in all L q ( 12) with q < oo. The extra re¬ 
sult of uniform Holder regularity of v(-,t) implies convergence in L°°(12) by standard 
interpolation results between L 1 and C a . □ 

Remark. We have only found a related result for the asymptotic behaviour of the 
present Dirichlet problem, recently announced in [41], that provides a first rough estimate. 
Indeed, their Theorem 2.6 says: 

• Let q > p. Let u(x,t) be the solution of the Dirichlet problem for the initial datum 
Uo G L°°(l l), if q > p and Uq G L 2 (12) if q = p. Then, the L q -norm of the solution goes to 
zero as t oo since we have the following estimate: 

IHOIILto) < (4.2) 

5 Positivity. Reflection principle 

The proof of Theorem 12.21 has two steps. In the first we prove the result for the case 
where 12 is a ball and Uq is radially symmetric. In the second we obtain the whole result. 

• Radial monotonicity. Before we prove positivity of the solutions when 12 is a ball is 
Uq(x) is radially symmetric around the center of the ball and decreasing with the distance, 
we need to prove a monotonicity principle in the radial direction. We may assume that 
12 = Br( 0). We know that in that case the solutions are radially symmetric. We want 
to prove that they are monotone in the radial direction. The comparison is done by an 
application of Aleksandrov’s reflection principle. This has been performed in the case of 
the fractional porous medium in j-571 Section 15], and we only have to repeat the process 
followed there. 
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A main reduction is that we only need to prove the reflection principle for the elliptic 
steps of the implicit time discretization, which amounts to prove the result for solutions 
of 

C P)S u + u = f. (5.1) 

posed in 12 with / E L 2 (12). We take a hyperplane H that divides into two half¬ 
spaces, and let 12i and H 2 the two pieces into which H divides 12. Let us assume that the 
symmetry II with respect to H that maps 12 1 into H 2 (i. e., H does not pass through the 
origin, which is contained in 12 2 ). The statement we need to prove is as follows. 

Theorem 5.1 Let u be the unique solution of (j5. 1 j) with data f E L 2 (12) ? f > 0. In the 
above circumstances and under the further assumption that 

f(x) < /(n(x)) in Oi (5.2) 


we have u(x ) < w(II(:r)) in Oi. 

The proof does not differ essentially from the one in [57]. As a corollary, we conclude 
that whenever / is radially symmetric and nonincreasing in |x|, so it u(x). Passing to a 
similar statement about the parabolic evolution is immediate. 

• I. We can now prove positivity for the solutions of the parabolic problem when 12 is a 
ball if u(x,t ) is radially symmetric around the center of the ball and decreasing with the 
distance (at all times). 

Suppose that u{r,t) is a radially symmetric solution, r = |x|, and it is nonincreasing in 
r. If for some t > 0 have u(ri,t) = 0 a point rq < R we will also have u(r,t\) = 0 for 
all r\ < r < R. By the time monotonicity we will also have u(r,t ) = 0 in the cylindrical 
region: {rq < r < R, 0 < t < fi}. Now, in that region u t — 0 while C P)S u < 0, since the 
points are points of minimum of u, and u is not identically zero. This is a contradiction 
with equation (11.11) . hence u cannot touch zero. We conclude that u(r,t ) > 0 everywhere 
in Q. 

• II. Positivity for general solutions. We take now a general domain and only assume on 
u 0 that is strictly positive in some ball, say u 0 (x) > s > 0 in B R (x 0 ) C LI. We consider 
the solution Ui(x,t) with domain B R (x 0 ) and initial data u 0 1 smooth, radially symmetric 
around xq, decreasing in r = \x — xq\. By the previous item, u±(x, t) is positive everywhere 
in Q(R-,x 0 ) := B R (x 0 ) x (0,oo). A simple application of the maximum principle implies 
that u(x,t ) > u\(x,t) in Q(R]x 0 ). Notice that maximum radius R we may take the 
distance from xq to the boundary <912. 

In order to propagate this positivity everywhere in Q = 12 x (0, 00 ), we select a point 
xi E 12 and a time T > 0 and we argue as follows: if |xi — xo| is less than the distance 
from xo to the boundary we are already done. If this is not the case there is a finite set of 
points yi, y 2 , ■ ■ ■ , y N = ^ 1 , such that \y k+1 - y k \ < \d k , where d k = dist (y k , <912). Picking 
times t 0 = 0 < ti < t 2 <■■■< t n = T and applying iteratively the previous argument, 
we conclude that u(x, t k ) > e k > 0 for x E B dk / 2 (y k ). In the end u{x 1 , T) > e N > 0. 
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Filially, when m 0 is not known to be positive in a ball, we wait a time e > 0 and apply 
the regularity results to show that u(x, e) must be positive in a ball since it is continuous, 
nonnegative and nontrivial. The rest of the argument is the same for t > e > 0. Since £ 
is arbitrary, we get positivity everywhere for all t > 0. Q 

6 Sharp boundary behaviour, asymptotics with rate 

In a remarkable paper m Ros-Oton and Serra studied the regularity up to the boundary 
of solutions to the Dirichlet problem for the fractional Laplacian, case p = 2. They con¬ 
sider the solutions u of the elliptic equation (—A ) s u = f with 0 < ws < 1 and / e L°°(fl), 
and satisfying Dirichlet condition u = 0 for x fL 0. They prove that u G C S (M. N ) and 
u(x)/d(x) s is C a in D and up to the boundary, for some a(0,1), where d(x) = dist(x ; dQ). 
They do that by developing a fractional analog of the Krylov boundary Harnack method. 
This assertion must be strengthened with the assertion that the quotient u(x)/d(x) s is 
strictly positive, which is done in Lemma 3.2 of the paper, see also j46j. The results were 
extended to the fractional heat equation in [33] . 

• We conjecture that this boundary behaviour result is true for p > 2 with the same power, 
u(x,t ) ~ d(x) s . The upper bound for F(x)/d(x) s is true by the already mentioned upper 
boundary behaviour estimate, [35], and as a consequence of the uniform bound in terms 
of U(x,t), the same upper bound applies to any solution u(x,t). But the sharp lower 
bound is missing at this moment, and it has important consequences. 

If this is accepted as a sensible hypothesis, we can prove the following version of the 
sharp asymptotic convergence with rate. 

Theorem 6.1 Let the conjecture hold and letu be nonnegative solution u(x,t) of problem 
DP with nontrivial data. Then there is a constant T > 0 that depends on the data such 
that 

U(x,t + T) < u(x,t) < U(x,t) (6.1) 

holds in D if t is large. This means that as t —> oo 

t 1 T p ~ 2 l u (x,t) = F(x) + 0(l/t). (6.2) 

Remark. The example u(x,t) = U(x,t + T) shows that the result is sharp. 

Proof. We only need to establish the lower bound, U(x, t+T) < u(x, t). By the maximum 
principle, this can be done at any given time, say t — t\ and then it will hold for all times 
t>t\. Now, at t — ti we have by the assumed regularity 

u(x,ti) > Cd(x) s , C = C(u,ti) 

while away from the boundary u(x,ti) is uniformly positive by continuity and the pos¬ 
itivity result of the previous section. We conclude that there is an £ > 0 such that 
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u(x,ti) > £ F(x). On the other hand, for T large enough 

U(x,t) = (h + T)- 1/(p - 2) F(ir) < eF{x) . 
The conclusion follows. □ 


7 Comments and open problems 


We list here a number of comments and questions that might interest the curious reader. 

• The difference between the action of operator C PtS in and on a bounded domain 
can be seen when we compute the lack of conservation of mass for a strong nonnegative 
solution of problem DP in a bounded Q. From the equation we have 


d 

dt 


u 


{.x , t ) dx 



£ PtS u(x, t) dx. 


After putting 1 q(o;) = 4>(x) and using the fact that u(x,t ) can be extended by zero if 
x ^ hi, the last integral equals, a 


r N 


jC PtS u(x,t)ladx = 


\u(x) - u(y)\ p 2 (u(x) - u(y)) (<f>(x) - <j>(y)) 
\x — y\ N+s P 


dxdy. 


We only need to take into account the case where i6ll and y ^ (and the converse), 
otherwise cft(x) — (p(y) = 0. Therefore, we get 


d 

dt 



u(x,t ) dx 



u(x)\ p 2 u(x) 
\x - y\ N+s P 


dxdy , 


which is never zero, and determines the loss of mass. In M. N such an integral does not 
exist and mass is conserved for all solutions for which the calculation is justified. 

• The regularity of strong solutions of Problem DP is a very ongoing interesting question. 
See the works [T6, 37] for non-degenerate problems of this type, where C a regularity of 
bounded solutions is proved. For fractional porous medium equations it has been proved 
in j2n> M] and |3H] for the different models. As we have noted above, sharp boundary 
regularity is really needed. 

• The universal upper bound is true for signed solutions, that are naturally produced by 

the theory when n 0 is not necessarily negative. We get in particular the following version 
for a signed solution: | u(x, i)| < This translates into corresponding bounds 

for the large-time decay of the form \u(x,t) \ = but the sharp versions of the 

asymptotic behaviour are lost since F(x) is not representative, generally speaking. See 
the analysis for the standard PME and PLE in pufj . 

• Theorem 13.21 states that the universal upper bound applies to all solutions with initial 
data in L 1 (h2). By passage to the limit it would also apply to solutions having as initial 
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data a nonnegative Radon measure. We will not enter into the theory for such solutions, 
since more basic questions are still open. 

• The universal upper bound and universal asymptotic behaviour are true for all p e 
(2, oo), but they do not hold for the limit values of p. Thus, they cannot be true for the 
limit case p = 2, since it is a linear equation and the size of the solutions must depend 
linearly on the size of the data. For the other limit, p —> oo, an evolution theory has to 
be developed, see [27], [39] for the proposed elliptic operator. The reader will notice that 
the time factor t -1 /h~ 2 ) of the Friendly Giant goes to 1 in the limit, its decay disappears. 

• Regarding the asymptotic behaviour in the linear case p — 2, it is known that solutions 
of the Dirichlet Problem decay exponentially, and more precisely, the large time behaviour 
of a solution takes the classical form u(x,t ) ~ c(u 0 )(p(x)e Xlt where 0(x) > 0 is the first 
eigenfunction and Ai > 0 is the first eigenvalue of £ 2 ,s = (—A) s . So the asymptotics is 
quite different. 

• The large time situation is even more different when 1 < p < 2. Though we will 
not study that range in this paper since it deserves proper attention, let us just point 
out that the method of Section [3] shows in a simple way that there a phenomenon of 
complete extinction in finite time, as the one already described for fast diffusion and fast 
p-Laplacian in the literature, [53]. A brief argument is as follows: we introduce a rescaled 
solution w(x, r) by means of new formulas 

v(x,t) = (T-t)- 1/{2 - p) u(x,t), t — —^— log(l/ (T — t)), (7.1) 

z p 

where T > 0 is arbitrary and the ranges of time are 0 < t < T and tq < t < +oo. This 
change leads to the equation for v: 

d T v + (2 -p)C p , s (v) = v. (7.2) 

Like in Section [3] we find that there is stationary positive supersolution W(x) for this 
equation. If moreover the data or the solution are bounded, i.e., if v(x,Ti ) < W(x) for 
some Ti, we conclude that v(x,t) < W(x) for all x G 12 and T\ < r < oo. Translating 
this for u we get 

u(x, t) < (T — t) l ^ 2 ~ p " > W{x) for xefl,ti<t<T, 

which implies that u(-,y) converges uniformly to zero as t /* T. We conclude that the 
orbit vanishes identically after a time T(u ) which is equal or less than the constant 

T of the supersolution. Analyzing the behaviour for t near the extinction time is an open 
problem. 

• The problem with right-hand side, d t u + £ PtS u = f is interesting. Existence and 
uniqueness of solutions is granted by the usual theory. See [33] for recent results on the 
elliptic problem. 

• Some authors have considered operators like £ PtS with more general nonlinearities than 
our <F(z) = \z\ p ~ 2 z. Our results do not apply to them, at least the sharp ones. 
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• There are alternative definitions of the fractional Laplacian on bounded domains in the 
linear case p = 2, that have been discussed recently in the literature, see e. g. [9] 021 150 j. 
No such alternative definition seems to apply in our p-Laplacian case. This is an intriguing 
question. 

• It would be nice to know more semi-explicit solutions of equation (II.ID . 

• Note that the evolution generated by Problem DP in fi, started with initial condition 
the eigenfunction (I3.7j) of j3U], produces in the rescaled u-formulation a natural connection 
between this profile and the friendly giant profile F(x). 

• The methods we have used do not allow us to settle the fine asymptotic behaviour of the 
solutions of the Cauchy problem posed on M. N . The idea is to construct the corresponding 
Barenblatt solutions, but such process is not easy. 

• Can we say that our p-s Laplacian equation with p > 2 is a degenerate parabolic 
equation when 0 < s < 1, as it is said of the same equation when s = 1? 
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